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Electromagnetic properties of periodic two-dimensional sub-wavelength structures consisting of 
closely-packed inclusions of materials with negative dielectric permittivity e in a dielectric host 
with positive eh can be engineered using the concept of multiple electrostatic resonances. Fully 
electromagnetic solutions of Maxwell's equations reveal multiple wave propagation bands, with the 
wavelengths much longer than the nanostructure period. It is shown that some of these bands are 
described using the quasi-static theory of the effective dielectric permittivity e qs , and are independent 
of the nanostructure period. Those bands exhibit multiple cutoffs and resonances which are found to 
be related to each other through a duality condition. An additional propagation band characterized 
by a negative magnetic permeability develops when a magnetic moment is induced in a given nano- 
particle by its neighbors. Imaging with sub-wavelength resolution in that band is demonstrated. 

PACS numbers: 41.20.Cv, 42.70.Qs, 42.25.BS, 71.45.Gm 

Electrostatic resonances of isolated nanoparticles have recently attracted substantial interest because of the in- 
t-^ , triguing possibility of obtaining very strong and localized electric fields. Applications of such fields include nanoscale 
biosensors pi Q , nanolithography 0, j and nonlinear spectroscopy 0, Q . Resonances occur when a metallic or di- 
electric particle with a negative frequency-dependent dielectric permeability e(u>) < is imbedded in a dielectric host 
Ci . (including vacuum) with a positive dielectric permeability eh > 0. The wavelengths A of the incident electromagnetic 
radiation that resonate with a small particle of a characteristic size d <C A depend on the particle shape and the 
functional dependence of e(u>). By changing the shape Q and internal composition Q of nanoparticles resonances 
can be shifted to the wavelength optimized for a particular application. Close proximity of other small particles 
can also strongly affect the resonances. For example, disordered conglomerates of nanoparticles have recently been 
q ■ analyzed using the quasi-static theory |8| and shown |9( to exhibit multiple coherent " hot spots" that can be used for 
transporting electromagnetic energy. Surface-enhanced Raman scattering (SERS) from molecules inside multi-particle 
aggregates has also been shown to be greatly enhanced 0, 0, . 

In this Letter we demonstrate how electrostatic resonances of periodic arrays of closely-spaced particles can be 
used for engineering propagation properties of electromagnetic radiation through such nanophotonic crystals. The 
emphasis of this work is on subwavelength photonic crystals (SPC) that have a periodicity d significantly smaller 
than the radiation wavelength in vacuum A = 2t: c/lu. By numerically solving Maxwell's equations, we identify two 
classes of waves supported by an SPC: (a) hybridized Dipole Modes (DM) that are characterized by a quasi-static 
period- independent dielectric permittivity e qs (oj), and (b) hybridized Higher-order Multipole Modes (HMM) that 
depend on the crystal period d. Two types of DM's are identified: almost dispersionless (non-propagating) collective 
plasmons (CPL) satisfying the uj{k) = uj^ dispersion relation (where w'f' are multiple zeros of e qs ), and propagating 
collective photons (CPH) satisfying the k 2 c 2 = u> 2 e qs (uj) dispersion relation. The mean-field dielectric permittivity e qs 
calculated from the quasi-static theory 0, 01 is found to be highly accurate in predicting wave propagation even for 
SPCs with the period as large as \/2ir. DM wave propagation bands are "sandwiched" between multiple resonance 
uJf^ and the cutoff frequencies of the SPC. For two-dimensional SPCs with a high lattice symmetry (square 
and hexagonal) a duality condition expressing a simple one-to-one correspondence between the resonant and cutoff 
O • frequencies is derived. 

The new HMM propagation bands are discovered inside the frequency intervals where e qs < and, by the mean-field 
• i-h , description, propagation is prohibited. HMM bands should not be confused with the usual high order Brillouin zones 
■ of a photonic crystal because the latter do not satisfy the d <C A condition. One HMM band defines the frequency 
range for which the sub-wavelength photonic crystal behaves as a double-negative metamaterial (DNM) that can be 
described by the negative effective permittivity e e fj < and permeability /i e // < |l3j| . Magnetic properties of 
the DNM are shown to result from the induced magnetic moment inside each nanoparticle by high-order multipole 
elecctrostatic resonances of its neighbors. It is shown that a thin slab of such DNM can be employed as sub- wavelength 
lens capable of resolving images of two slits separated by a distance <C A. 

To start, consider a TM-polarized electromagnetic wave, with non- vanishing H z , E x , and E y components, incident 
on an isolated dielectric rod (infinitely long in the z— direction) with e(u>) < 0. The incident em wave is strongly 
scattered by the rod when its frequency lo coincides with that of the surface plasmon found by solving the nonlinear 
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FIG. 1: Two strongest hybridized electrostatic dipole resonances of the 2-D square lattice of plasmonic rods contributing to 
e qs : lu^ = 0.38^ (left); u>2 = 0.63oj p (right). Shown are the potential isocontours and the electric field arrows. Spatial 
dimensions are measured in 2c/u) p units. SPC parameters: d — c/uj p and R = 0.45d. 



eigenvalue equation for H z : 



^H Z )=^H Z1 (1) 



where H z — > far from the rod. The electric field of the surface wave is given by E(x) — — i[c/u>e(x)]e z x Vi? z , 
where e(x) — eh outside and e(x) = e{uS) inside the rod. In what follows we assume that the rods are in vacuum, 
i. e. et = 1. Note that the lines H z = const are the electric field lines. For a sub-wavelength rod the rhs of Eq. JIJ can 
be neglected. Moreover, \E\ ^> \H\, and the description using the electrostatic potential <j> is appropriate: E = —Vcj>. 
Hence, two equations are simultaneously satisfied for a sub-wavelength rod: 



-V • [e- L VH z j =0 and - V • [eV^j = 0. (2) 

Equations (J2J illustrate that in electrostatics there are two equivalent descriptions of the electric field: potential 
description and field line description. Note from Eqs. © that if a surface wave is supported by a rod of an arbitrary 
transverse shape for e = e\, then a surface wave is also supported for e = 1/ei Q. 

For a round cylinder of radius R the surface-wave multipole solutions of the second of Eqs. J2J are given by 
^out = ( r / R)~ m e x P (imO) outside and <j)\™^ = (r / R) m exp (imff) inside the cylinder, and m > I. The continuity of 
ed<f>/dr is satisfied for every m > 1 for cjW such that e(uW) = —I. This de generacy of the multipole resonances is 
specific to the 2-D lattice of round cylinders. This property makes the square SPC of round cylinders particularly 
amenable to engineering its electromagnetic properties using electrostatic resonances: various multipoles of a given 
rod are strongly hybridized by the proximity of other rods. This gives rise to a rich and easily controllable set of 
hybridized electrostatic resonances whose frequencies depend on the areal fraction occupied by the cylinders. Earlier 
approximate calculations 01 of the three-dimensional lattice of metallic spheres failed to demonstrate the type of rich 
multi-resonant band structure found for 2-D structures and shown in Fig.[21for two reasons: (a) there is no frequency 
degeneracy in three dimensions, and (b) the volume fraction occupied by metallic spheres is 3-D is smaller than in 
2-D. 

Note that there is no monopole (m = 0) electrostatic resonance (although there is an electromagnetic Mie resonance 
for rods with a very high positive e |15|1 for an isolated cylinder. Such a resonance would reveal that the corresponding 
solution for H z has an associated azimuthal ^-independent current around the cylinder and, thus, a non-vanishing 
magnetic moment M — (l/2c)(r x J), where the average is taken over the unit cell. However, as shown below, the 
octupole (m = 4) electrostatic resonances in a square lattice of closely-packed cylinders hybridize in a way of inducing 
a resonantly excited magnetic moment. This magnetic moment manifests itself as a negative effective permeability 
fJ-eff of the structure, and an additional propagation band of DNW's in a narrow frequency range. 

In the rest of the paper we concentrate on a specific SPC: a square lattice of round (R = 0.45(f) plasmonic cylinders 
with e(ui) = 1 — Up/oj 2 characteristic of collisions-free electron gas, lattice period d = c/u) p . Very similar results are 
expected for polaritonic rods with e(u>) — eooiuj 2 — oj\q) / (u> 2 — luxo)^ w ith e < for u>to < lu < ^lo- To this SPC we 
apply the standard procedure for calculating the quasi-static dielectric permittivity e qs (uj), and later compare 

the band structure described by e qs (u>) to that obtained by solving the fully-electromagnetic Eq. ^ 

The material-independent e qs is calculated 0, as 

N Fi 

e qs = l-py2 ( 3 ) 
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FIG. 2: (a) Frequency dependence of the quasi-static dielectric permittivity e q3 from Eq. 10. Vertical lines: resonances, 
horizontal axis intercepts: cutoffs, (b) Collective photons (CPH) supported by the SPC: theoretical prediction of the quasi- 
static theory, k — ^/e q3 Lj/c (solid lines) and numerical simulation of Eq. (0 (circles and diamond). Note the backward- wave 
mode (dio/dk < 0) marked by diamonds which is not described by the quasi-static theory. SPC parameters as in Fig.Q 



where, for the plasmonic rods in vacuum, s(u>) = [1 — e(u>)] 1 = to 2 /ui 2 , s, = [1 — ei\ 1 is the i'th out of N > 1 
hybridized dipole resonances, Fi its oscillator strength calculated below, and p = tt R 2 /cP. Electrostatic resonances Sj 
are found by solving the eigenvalue equation for the potential eigenfunctions fa inside a unit cell of the structure: 

(9(f) V^j = stf 2 fa, (4) 

where fa satisfies the following homogeneous conditions at the unit boundaries (x,y) = (±d/2, ±rf/2): (a) fa and its 
derivatives are periodic; (b) fa(x = ±d/2) = 0; (c) d y fa(y — ±d/2) — 0. Equation (fJJ follows from Eq. (J2J). Physically, 
these eigenfunctions describe the electric potential distribution when a vanishing ac voltage (with frequency u> such 
that e(w) = £i) is applied between x — ±d/2 capacitor plates. The capacitance of such an imaginary capacitor, 
equal to the ratio of the charge to the voltage drop, is given by C — e qs d, and becomes infinite according to Eq. (J3J - 
Another eigenfunction fa corresponding to the voltage drop between y — ±d/2 plates is obtained by a 90— degree 
spatial rotation of fa. 

Because the square lattice is invariant with respect to the transformations of the Ca v point group [Tfil | , all periodic 
solutions transform according to one of the irreducible representations (irreps) of C^ v : four singlets (commonly labeled 
as Ax, A 2 , Bi, and B 2 ) and one doublet E. The electrostatic eigenfunctions 4>i and fa have the symmetry of E. 
Inside a given rod each fa can be expanded as the sum of multipoles: fa{r,9) = J2Zo A\ 2l+1) (r / R) 2l+1 cos (21 + 1)9. 
A straightforward calculation following Ref. yields the oscillator strength proportional to the dipole component 

of fa: Fi — {A\ ) 2 / Y^aLq{^ + 1) ■ For our structure there are three significantly strong hybridized dipole 

(^-symmetric) resonances : (si = 0.1433, F x = 0.8909), (s 2 = 0.4025, F 2 = 0.064), and (s 3 = 0.6275, F 3 = 0.0366). 

The plots of the two lowest resonances are shown in Fig. The first resonance is primarily dipolar (ex cos fa) while 
the second one has a significant sextupolar (oc cos3</>) component. Thus, the close proximity of the rods in the lattice 
results in a strong hybridization of the odd multipoles with the dipole. Moreover, the hybridized dipole resonances 
uj^ = 0.38w p , — 0.63w p , and u>^ = 0.79w p occur at the frequencies controllably different (through the R/d ratio) 
from that of an isolated rod, lo^ = w p j\[2. 

The frequency dependence of the e qs for a plasmonic material is plotted in Fig. [2a). The vertical lines in Fig.|2fleft) 

are the plasmonic structure resonances which occur at oj^ — 0.38u p , uj^ = 0.63w p , and oj^ = 0.79o-> p . Quasi-static 

theory predicts no propagation in the region of e qs < 0. Three frequencies oj\ c ^ for which e qa {u)\ ) — are called 
the cutoff frequencies. From Fig.^a), there are four propagation bands (where e qs > 0) allowed by the quasi-static 
theory for u) < uj p . The quasi-static theory predicts that this plasmonic SPC acts as an effective medium supporting 
two types of waves: four collective photons (CPHs) satisfying the \k\ 2 = e qs ui 2 /c 2 dispersion relation (two of which are 

indicated by a solid line in Fig.^b)), and three non-propagating collective plasmons (CPLs) satisfying the uj(k) — uj^ 
dispersion relation (not shown) . Numerical solutions of the fully electromagnetic Eq. (JTJ described below confirm these 
conclusions, yet reveal additional modes whose behavior is not described by the quasi-static e qs . 

An important duality principle for the cutoff and resonance frequencies exists for square and hexagonal lattices 
(correspondingly invariant with respect to the C^ v and Cq v point group transformations): for each resonance frequency 

uj^-p there exists a cutoff frequency lj^ for which s f w | C ^ = 1- This principle follows from an observation that, at 



4 



the cutoff, H z satisfying the first Eq. obeys the same boundary conditions as <f> and <j>. Hence, from Eq. J5J follows 
that for each resonance frequency uj^ there exists a cutoff frequency to^ for which e(oj^) — l/e((jj). Numerical 
calculation of the zeros of e(s) from Eq. © indeed confirms the duality principle. 

To verify that e qs is sufficient for accurate description of em wave propagation through the plasmonic SPC, the fully 
electromagnetic Eq. was numerically solved as a nonlinear eigenvalue equation for ui 2 /u>p for different wavenumbers 

k = e x k, and the resulting dispersion relation u>(k) plotted in Fig.^b) for three propagation bands. Although d — 1 /uj p 
is not infinitesimally small compared to the radiation wavelength, it is apparent from Fig. [21 that the numerically 
calculated points (circles) accurately fall on the solid lines predicted for the CPH's by the scale- independent e qs . The 
essentially flat propagation band oj(k) 0.61o-> p not shown in Fig. H£b) also agrees with the CPL dispersion relation 
obtained from e qs . Therefore, for several frequency bands, the plasmonic SPC indeed is an effective medium described 
by the scale- independent e qs . 

Importantly, Fig. (2fb) reveals that there is another propagation band (diamonds) in the frequency range for 
which no propagation is expected due to e qs < 0. Note that the mode's group velocity v g — duj/dk < op- 
poses its phase velocity - an indication that we're dealing with a DNM. For k = this mode's H z has the 
symmetry of the A\ irrep of the symmetry group Cav, and can be expanded inside a given plasmonic rod as 
H z (r,9) = J2T=o^ 4k ^ [l<ik(wV— er/^/hkiwy/^eR/c)] cos (4fc0), where Ii is the modified Bessel function of order 
I. Because there is no dipole component in H Zl the A\ mode does not contribute to the quasi-static permittivity 
e qs . For the SPC at hand, the largest term in the expansion is the octupole term A^', and the next largest is the 
monopole term A^ that is responsible for the magnetic moment induced in the photonic structure as explained 
earlier. Therefore, the mode is an HMM, with predominantly m = 4 component. Because of the finite magnetic 
moment, a single quantity e qs cannot describe the HMR resonances, and two frequency-dependent parameters are nu- 
merically evaluated: dielectric permittivity e e // and magnetic permeability [i e ff- The procedure for expressing these 
effective quantities for a periodic structure using the cell-averaged electric and magnetic fields has been described 
elsewhere [TH Il7| . 

For two dimensions, and assuming that the elementary cell of the SPC is centered at the origin, we intro- 
duce several variables by averaging H z and E over the sides or the area of the elementary cell of the photonic 
crystal: B z = dr 2 j dAH z (x,y), H z = H z {x = -d/2,y = -d/2), E y = d' 1 dy E y {x = -d/2,y), and 

D y = d^ 1 /_^/2 dxE y (x, y = —d/2). For k = ke x Maxwell's equations in the integral form become kcg(kd)E y = ujB z 
and kcg{kd)H z — ujD y , in complete correspondence with Maxwell's equations in the medium for which B z = (i e f fH z 
and D y = e e ffE y . Dimensionless factor g(x) = i[l — exp (ix)]/x —> 1 for kd -C 1 is the slight modification accounting 
for non-vanishing lattice period. 

The magnetic permeability /x e // is affected because the mode carries the electric current which produces a finite 
magnetic moment. The magnetic nature of the A\ mode is due to the non-vanishing coefficient of the monopole 
A (o) term in the multipole expansion. The monopole is responsible for the ^-independent component of the az- 
imuthal electric field Eg = —i[c/uje(r)]d r H z . The corresponding electric current in the negative-e rod given by 
Jg = — A^uj/c(1 — l/e)Ii(ujr/c)/Io(ujr/c) produces a magnetic moment density M = (l/2c)(f x egjg), where the 
average is taken over the unit cell, and can be shown to be M = — (pAo/An)(l — l/e)l2(u>R/c)/Io(u)R/c). 

The effective permittivity and permeability have been calculated for a range of wavenumbers k = ke x and the 
corresponding frequencies u>(k). For kd <C 7r it follows from the analyticity of u>(k) that the frequency depends only on 
|fc| and not on its direction. For fco = 0.6/d and luq = Q.Qu p (or n e f / = —1) we numerically computed that /i e // = —2.35 
and e e ff — —0.427. Therefore, at this frequency our SPC is a DNM. Note that the hybridized monopole/octupole 
resonance affects not only the magnetic permeability of the SPC, but also the dielectric permittivity: the mean-field 
calculation using Eq. @ yields e ?s (wo) = —0.65 that is significantly different from e e //- 

DNM-based flat super- lenses capable of sub- wavelength imaging have been proposed ^| . The condition for super- 
lensing is that the DNM with the dielectric permittivity e < is embedded in a host medium with = — e. We 
have tested a six-period thick plasmonic SPC for the super-lensing effect by embedding it inside the hypothetic host 
with eh = 0.55. This particular choice of — e e // < th < — £<js was not optimized, and is one among the several that 
showed super-lensing. To verify the sub-wavelength resolution, we simulated the distribution of the magnetic field 
\E\ behind a screen with narrow slits of width A y = A/5 separated by a distance 2A y . As depicted in Fig. [31 where 
only two slits are shown, a planar wave with frequency uj = 0.6o> p is incident on the screen from the left. A six-period 
long plasmonic SPC of width D — 0.6A is positioned between < x < D. The distribution of \E\ in the x — y 
plane is shown in Fig.^a). Also, in Fig.[3{b) \E\ is plotted in two cross-sections: the object plane right behind the 
screen (at x = —D/2 + A/10, solid line), and in the image plane (at x = 3D/2 — A/10, dashed line). Object plane is 
slightly displaced from the screen to avoid £?-field spikes at the slit edges. The two sub-wavelength slits are clearly 
resolved. Increasing the incident frequency by just one percent (outside of the DNM band) results in the complete 
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FIG. 3: (a) Magnetic field distribution behind an illuminated periodic slit array, with a six-period SPC, parameters as in 
Fig- Q (b) \E\ in the object plane (solid line); in the image plane for ujo — 0.6lo p , without damping (dashed line) and with 
damping characteristic of silver (dot-dashed line); in the image plane for ui = 0.606w p (dotted line). 



loss of resolution in the image plane (dotted line). 

While the DNM band for the plasmonic SPC is quite narrow, [uj(k = 0) — u)(k = Tr/d)]/u>(k = 0) = 0.055, it is still 
broader than the collisional linewidth for some plasmonic materials. For example, for silver e = — uj^/uj(uj + ij), 
where €b ~ 5, uj p — 9.1 eV, and 7 = 0.02 eV 0|. Fig. El^b) (dash-dotted line) confirms that, although finite damping 
j/uip — 0.002 reduces the field amplitude in the image plane, it does not affect the image contrast. The band flatness 
in a plasmonic SPC translates into very sharp excitation resonances, and large enhancements of the incident field. 
For example, the very close proximity of two flat propagation bands in Fig. Bib) (second DM and the A\ HMM) can 
be exploited for maximizing the structure response at the incident and Raman-shifted re-emitted frequencies, which 
is essential for SERS. 

This work is supported by the NSF's Nanoscale Exploratory Research Contract No. DMF0304660. 
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